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1. Introduction 

Algebraic computing in relativity and gravitation dates back more than 
thirty years, but only relatively recently has hardware of sufficient power to 
tackle large scale calculations become commonplace. Whereas it is generally 
understood throughout the relativity community that there are a number 
of packages available, the diversity of problems that the available packages 
can help with is not so widely appreciated. This session was devoted to 
computer algebra for relativity and gravitation from the point of view of 
developers. In this summary I expand this to include some background and 
outline what is available for users in the field. 

2. Background 

If indeed "the beginning is the most important part of the work" 1 let me 
start by emphasizing the fact that modern computer algebra systems are 
easy to use and can do rather more than calculate the components of the 
Riemann tensor! About fifteen years ago it would have been reasonable to 
list active users of computer algebra systems in the relativity community 
(e.g., (Cohen et al, 1984)). This is no longer feasible. And yet, for example, 
we find " Algebraic Computing in General Relativity" but one of six topics 
in "Numerical Relativity" within the newly established Living Reviews In 
Relativity. 2 To this day computer algebra is, I believe, an underused tool. 
It can certainly be argued that it is also under supported (Petti, 1994), and 
perhaps even misunderstood. Three questions come to the fore: How did 
computer algebra develop, how is it distinct from more traditional uses of 



ipiato, The Republic. Book II. 377B 
2 The web address is www.livingreviews.org 
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the computer, and what can it do for researchers in relativity and gravita- 
tion? 

As regards the history, useful references are (D'Inverno, 1980) and (Mac- 
Callum, 1996). In computer algebra calculations, expressions tend to inflate 
before they simplify. 3 This intermediate expression growth dictates garbage 
collection, a standard feature of LISP and C but not of FORTRAN. More- 
over, since the final size of expressions are not known a priori, recursive 
algorithms, which are impossible in FORTRAN, are an essential feature 
of any computer algebra system. 4 The consequence is that at present the 
available algebraic engines are LISP-like (e.g., REDUCE 5 and Macsyma 
6 ) or C-like (e.g., MapleV 7 and Mathematica 8 ). Intermediate expressions 
must of course be stored and simplified. The key requirements of a com- 
puter algebra system are; in the software the ability to simplify, and from 
the hardware ample high speed memory (RAM). CPU speed is of some 
interest, but is considerably less important than it used to be. 9 

The environment that a computer algebra system presents to a user, fa- 
miliar with numerical computing, is strange territory, though perhaps not 
at first sight. The level of complexity of a problem submitted to the system 
is not always clear beforehand. Consider, for example, a calculation of the 
number of terms in the covariant Ricci tensor in n dimensions for a general 
form of the metric (that is, n(n+l)/2 functions of n variables) (Fee et al, 
1988). With n=4 the number of terms in each diagonal component of Ricci 
is 9,990 and in each off-diagonal component there are 13,280 terms. This 
once classic problem can now be handled on almost any PC with appropri- 
ate software. However, the number of terms grows rapidly with n. At n=5 
the number of terms in each diagonal component of Ricci rises to 298,134, 
and to 410,973 in each off-diagonal component. The case n=5, without re- 
course to custom code, requires a dedicated resource of about 1/2GB of 
RAM. 10 A increase in difficulty from n=4 to n=5 is certainly anticipated, 
but the magnitude of this change is difficult to predict. Another feature es- 

3 An example of this intermediate growth is provided by the Tomimatsu-Sato space- 
times. A discussion can be found at www. astro. queensu.ca/'lake/demo. html and at 
www.scar.utoronto.ca/~harper/redten/node60.html 

4 See e.g., (MacCallum, 1996). 

5 www.rrz. uni-koeln.de/REDUCE 

6 www.macsyma.com 

7 www .maplesoft.com 

8 www.wri.com 

9 A reasonable number to gauge hardware speed for computer algebra is SPECint95 
(www.specbench.org/osg/cpu95/results/cint95.html). At the time of writing, PCs are 
comparable to traditional workstations on this measure and the associated SPECint95 
numbers have increased by about a factor of four in the last year and a half. 

10 At the present time, within MapleV, this calculation also requires a 64-bit architec- 
ture, a feature one can expect of PCs in the near future. 



COMPUTER METHODS IN GR: ALGEBRAIC COMPUTING 



3 



sentially universal to all modern computer algebra systems is the need for 
effective use of the available simplification procedures. This requires some 
practice, though some general guidelines are usually available for particu- 
lar systems. 11 For most situations these general rules will give adequate 
performance 12 , and reduce the calculations for even complex spacetimes 
to an essentially trivial exercise, though one that could not be attempted 
by hand. However, failure to follow the rules can make even simple calcu- 
lations intractable. In practice, computer algebra calculations tend to be 
either trivial, or impossible. A computer algebra system is a tool, and the 
skill of the user is measured by the ability to turn the impossible into the 
trivial. In reality this skill is rather easy to obtain. 

What is available, and what features can a user hope to find? In an 
attempt to answer the first question I maintain a list. 13 The second question 
is harder to deal with. I outline below some thoughts on what one would 
like to see in a modern package: 14 

— A transparent easy to use feature-rich interface; 

— A full slate of simplification and manipulation routines; 

— The ability to work in a variety of formalisms with seamless integration; 

— A menu which contains all standard objects along with derivatives in 
any index configuration, operators, support for specialized symmetries, 
and Petrov and Segre typing; 

— The ability to define, in a transparent way, new objects taking into 
account their symmetries; 

— The ability to work with multiple metrics simultaneously and deter- 
mine junction conditions on timelike, spacelike and null boundaries; 

— The ability to develop customized packages (in any number of dimen- 
sions) that can be used in conjunction standard functions and objects; 

— An on board database of known solutions of Einstein's equations with 
"equivalence problem" tools; 

— A robust tightly integrated abstract package; 

Unfortunately, no single package is currently available which does all of 
the above. 



11 For example, www.astro.queensu.ca/~lake/simp.html 

12 Usually, it is the answer that is of interest and not the fact that the simplification 
strategy is optimal. When optimal simplification strategies are the prime concern the 
problem is more involved because of the large number of simplification procedures that 
are usually available and the size of the resultant parameter space to be explored. 

13 www.astro.queensu.ca/~lake/course/algebra.html 

14 A personal note of caution is in order here. The concept that you get what you pay 
for does not seem, at least to me, to extend to the software "industry". 
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3. Presentations 

At the session the "equivalence problem" (the unique coordinate-free char- 
acterization of a space, see e.g., (MacCallum and Skea, 1994)) once again 
dominated the contributions. Jan Aman (Aman et al, 1997) discussed a 
class of Riemann-Cartan Godel-type space-times from the point of view 
of equivalence problem techniques by using a suite of computer algebra 
programs called TCLASSI. A coordinate-invariant description of the gravi- 
tational field for this class of space-times was presented and the possible 
groups G r of affine-isometric motions were discussed. Conditions for space- 
time (ST) local homogeneity were derived and the local equivalence of the 
ST homogeneous Riemann-Cartan Godel-type space-times was discussed. It 
was shown that they admit a five-dimensional group of affine-isometries and 
are characterized by three essential parameters £,m 2 ,uj: identical triads 
(£, m 2 ,uj) correspond to locally equivalent manifolds. The algebraic types of 
the irreducible parts of the curvature and torsion tensors are also presented. 
Denis Pollney (Pollney et al, 1997) described a suite of programs to help in 
the calculation of an invariant classification of a space-time in MapleV. The 
algorithm used is largely based on (though not equivalent to) the algorithm 
used in JanAman's CLASSI extension of the SHEEP computer algebra sys- 
tem. He described the implementation of the classification algorithms in 
MapleV within the GRTensorll package where factorization, integration or 
graphical facilities are available. The basic algorithms employed have been 
documented and are to be placed in the public domain. The method used 
for classification is that of (Aman and Karlhede, 1981) which uses algebraic 
properties of the curvature spinors to align the basis vectors and restrict 
the SL(2,C) freedom of the tetrad as far as possible in order to reduce the 
number of independent components to be calculated. As he pointed out, 
an important consideration in implementations of this method is the choice 
of invariant canonical forms for the curvature spinor components and their 
derivatives. A specific set of canonical forms for symmetric spinors was 
proposed, with a discussion of alternatives, with particular attention to 
computability, practical use, and interpretation. To this end, he described 
new tools in GRTensorll for the calculation and manipulation of symmetric 
spinors, including the calculation of their independent derivatives to any 
order as well as the ability to carry out general SL(2,C) rotations of their 
components. Ray d'Inverno (d'Inverno et al, 1997) presented an interac- 
tive database of exact solutions of Einstein's field equations available on the 
web. 15 The intention of the database is not simply to provide a list of exact 
solutions, but also to allow the user access to the invariant classification of 

15 Available at edradour.symbcomp.uerj.br, at www.maths.soton.ac.uk/~rdi/database 
and at www.astro.queensu.ca/~jimsk 
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these solutions. The database provides access to a unique characterization 
of the Riemann tensor and its covariant derivatives for over 200 metrics, 
which the user can view without performing calculations. The database 
allows three modes of search: (i) by one or more discrete characteristics 
which the user selects by radio buttons, (ii) by keyword search, and (hi) by 
equation number as listed in (Kramer et al., 1980). Future plans include an 
interactive interface to invariant classification routines to enable a user to 
generate invariant classifications "on-line". 

In addition to discussions of the "equivalence problem", Jim O'Connor 
(O'Connor, 1997) described use of the package Dimsym 16 (a program pri- 
marily for the determination of symmetries of differential equations) to 
seek collineations for the Kimura metrics. He demonstrate that when a 
Killing tensor is known for a metric one can seek an associated collineation 
by solving first order equations that give the Killing tensor in terms of 
the collineation rather than the second order determining equations for 
collineations. Anders Hoglund (Hoglund, 1997) described Tensign, a new 
program for indicial calculations of tensors. This program takes advantage 
of a cursor that you move around in the expression so that you can perform 
your operations just where you want them. One example of a calculation 
that has been done with Tensign is the wave equation for the Lanczos po- 
tential (Edgar et al, 1997). Another example is the investigation of the 
existence of the Lanczos potential in dimensions larger than four (Ander- 
sson et al, 1997). I described some innovative work by Peter Musgrave 
(Musgrave, 1997) which is in effect the first (as far as I am aware) vir- 
tual algebraic calculator for the web. GRLite 17 is an experimental easy 
to use GUI to GRTensorll. The Java applet 18 is a client which connects 
to a Java-based server. This server in turn starts a MapleV session and 
acts as an intermediary translating the button clicks and menu selections 
into GRTensorll commands. GRLite provides a prototype for a platform 
(hardware and algebraic engine) independent unified GUI to a variety of 
computer algebra programs on and off the Web. 

4. Some Recent Examples 

Some recent examples, selected from preprints received in the last two 
weeks, demonstrate how computer algebra has been used very recently as a 
tool. Unruh 19 has reexamined cosmological long wavelength perturbations 

16 www.latrobe.edu.au/www/mathstats/Maths/Dimsym/ 

1 7 www. astro . queensu . ca/ ~ lake / vac .html 

18 In Pune I was able to demonstrate GRLite (with the server in Canada) on a modest 
computer (Intel 386@33MHz) by way of the Kretschmann scalar for the Kerr-Newman 
metric. When a connection was available, the demonstration took but a few seconds. 

19 astro-ph/9802323, start at xxx.lanl.gov 
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and presented an exact solution to the long wavelength perturbations for 
the scalar modes and for a scalar field theory with arbitrary potential. The 
algebraic manipulations necessary in driving the equations and in reducing 
the Hamiltonian action were done with the help of a computer algebra sys- 
tem. The same author has also recently examined radiation reaction fields 
for an accelerated dipole for scalar and electromagnetic radiation. 20 The 
radiation reaction fields at the location of a point source for a massless 
scalar field, and the reaction fields at the location of a point dipole source 
for the electromagnetic field were derived using the Penrose integral. The 
extensive calculations were aided in an essential way with a computer alge- 
bra system. Glass and Krisch 21 have extend the Vaidya radiating metric to 
include both a radiation field and a string fluid. Assuming diffusive trans- 
port for the string fluid, they find new analytic solutions of Einstein's field 
equations. These represent an extention of Xanthopoulos superposition. To 
end, I would like to point out some recent work by George Davies. 22 He 
has developed computer algebra methods for calculating black hole pertur- 
bations which allow the "second order Zerilli function" 23 to be calculated 
in but a couple of moments on a PC. Work is in progress to extend this 
analysis to the Kerr metric. 
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